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Due to the ubiquitous presence of a Berry phase, in most cases of dynamical Jahn- Teller systems the 
symmetry of the vibronic ground state is the same as that of the original degenerate electronic state. 
As a single exception, the linear H ® h icosahedral model, relevant to the physics of Ceo cations, is 
determined by an additional free parameter, which can be continuously tuned to eliminate the Berry 
phase from the low-energy closed paths: accordingly, the ground state changes to a totally-symmetric 
nondegenerate state. 
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The traditional field of degenerate electron-lattice in- 
teractions (Jahn- Teller effect) in molecules and impurity 
centers in solids [[ljj|] has drawn interest in recent years, 
excited by the discovery of new systems calling for a revi- 
sion of a number of commonly accepted beliefs. Several 
molecular systems including Cgo ions, higher fullerenes 
and Si clusters, derive their behavior from the large 
(up to fivefold) degeneracy of electronic and vibrational 
states due to the rich structure of the icosahedral sym- 
metry group. H Novel Jahn- Teller (JT) systems have 
therefore been considered theoretically, disclosing 
intriguing features, often related to a Berry phase 

JTo| in the electron-phonon coupled dynamics. 

As it is well known, the molecular symmetry, reduced 
by the JT distortion with the splitting of the electronic- 
state degeneracy, is restored in the dynamical Jahn- Teller 
(DJT) effect, where tunneling among equivalent distor- 
tions is considered. The vibronic states are therefore la- 
belled as representations of the original point group of 
the undistorted system. In the weak-coupling regime, 
for continuity, the ground state (GS), in particular, re- 
tains the same degenerate representation as that labelling 
the electronic level prior to coupling. A priori, there is 
no particular reason for this to continue at larger cou- 
plings. However, it appears empirically j^j that in all 
linear DJT systems studied before the late nineties, the 
GS symmetry remains the same at all couplings. The 
explanation of this observation was a great outcome of 
the Berry-phase [[To) scenario: the phase entanglement in 
the electron-phonon Born-Oppenhcimer (BO) dynamics, 
H[Il-13 originating at electronically-degenerate high- 
symmetry points, seemed a universal feature of the DJT 
systems. 

In this context, it came as a surprise the discovery 
of the first linear JT system showing a nondegenerate 
GS in the strong-coupling limit: [^9) the spherical model 
2?( 2 ) ® S 2 ^ of electrons of angular momentum L = 2 
interacting with vibrations also belonging to an / = 2 
representation. This system turns out to be a special 
case of the H (g) h icosahedral model, for a 5-fold de- 
generate H electronic state interacting linearly with a 



distortion mode of the same symmetry h. |^],[l4| In that 
special case, it was shown that, for increasing coupling, a 
nondegenerate A excited state in the vibronic spectrum 
moves down, to cross the H GS at some finite value of 
the coupling parameter, thus becoming the GS at strong 
coupling. [^|J9| This phenomenon is a manifestation of the 
absence of Berry-phase entanglement in the coupled dy- 
namics. P,^5[ 

In this Letter we study the linear H ®h model in its 
generality. We analyse in detail the connection between 
the symmetry/degeneracy of the vibronic GS and the 
presence/absence of a Berry phase in the coupled dynam- 
ics. This model owns its peculiarities to the non-simple 
reducibility of the icosahedral symmetry group. In partic- 
ular, the H representation appears twice in the symmet- 
ric part of the Kronecker product of the H representation 
with itself: 



{H®H} {S) = a®g®h [1] ®h [2] 



(1) 



There are, therefore, two independent sets of Clebsch- 
Gordan (CG) coefficients 
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(H, mi; H, 777,2 1 h, m) 
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for the coupling of an H electronic state with an h vibra- 
tional mode, identified by a multiplicity index r = 1,2. 
p6| Of course, since the two h states are totally equiva- 
lent and indistinguishable, symmetry-wise, the choice of 
these orthogonal sets of coefficients has some degree of 
arbitrariness: the free parameter a in the combination 



(a) = cos a C, 



mi ,m2 



(3) 



The coefficient C™ x m (a) coincides with 
1 for a 



accounts for it 

the r = 1 and r = 2 values JT 

respectively. Also, for a — ~ arctan (3/\/5j = —a, 
becomes equivalent to the spherical CG coefficient. 

The basic Hamiltonian for the H ® h model can be 
written: 



= and a = ^ 
it 



H = H hlirm (huj) + H c - V (ghuj, a) 



(4) 
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At the boundary angles 



H c - V (ghu>,a) 
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(5) 



^ ] c t , mcl ni c- m2 C mi m2 (a) , (6) 
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where q m is the distortion coordinate (with conjugate 
momentum p m ) and cj^ is the electronic operator in stan- 
dard second-quantized notation. 

The novelty introduced by the a-dependent CG coeffi- 
cients reflects the fact that the group does not determine 
completely the form of the linear coupling as, for exam- 
ple, in cubic symmetry. The specific value of this an- 
gle must be established case by case by detailed analysis 
of the phonon mode and its coupling with that specific 
electronic state. Indeed, in a realistic case such as, for 
example, C^" ions, each h mode is characterized not only 
by its own frequency u)i and scalar coupling but also 
by its particular angle of mixing a,. 

For intermediate to strong coupling, the interesting 
nonperturbative regime, the customary framework is 
the BO separation of vibrational and electronic motion: 
when the splitting among the five potential sheets (pro- 
portional to g 2 ) is large, the electronic state can be safely 
assumed to follow adiabatically the lowest BO potential 
sheet, while virtual inter-sheet electronic excitations may 
be treated as a small correction. The BO dynamics is de- 
termined by the lowest eigenvalue of the interaction ma- 
trix 5 = ^2 QmV^ 71 ' in the electronic space. This matrix 
is obtained from (^J) by the same technique described in 
Ref. |lfj : since it is a simple generalization of that ob- 
tained for for 2?( 2 ) ® S- 2 \ |L5| here for brevity we report 
only the expression of the diagonal matrix elements of 
the matrix, corresponding to the coupling to a pure 
go distortion: 
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This form makes it clear that a shift a — > a + tt intro- 
duces a sign change in the coupling matrix, and it can 
be compensated by a reflection q — > — q. We will re- 
strict therefore, without loss of generality, to the interval 
< a < tt. 

The electronic eigenvalue ^1 cos a is the lowest for 
a < a s and a > tt — a s (region a): in this range the 
BO potential presents six absolute minima, one of which 
is lying along the qo pentagonal axis, with energy low- 
ering E'das = — g 2 /10 cos 2 a (in units of TuS). However, 
influenced by the y( m ^°) matrices, in the complemen- 
tary interval a s < a < tt — a s (region b), ten trigo- 
nal distortions become the absolute minima, with energy 



gain E clas 

(a = a s and tt — a 8 ), all pentagonal and trigonal minima 
become degenerate, and part of a continuous degenerate 
4-dimensional (4-D) trough || of depth i? c i as = -g 2 /28. 

We come now to the role of the Berry phase in this sys- 
tem. As well known, the geometrical phase is related to 
conical degeneracies of the two lowest BO potential sur- 
faces. [§ In the 2? (2) <g> system ||| (the a = tt - a s 
case of the model studied here) the flat minimum trough 
presents tangentially degenerate points. For that case, it 
was shown that the tangential contacts provide a mecha- 
nism for getting rid of the Berry phase. |l5|,[l8| For generic 
a instead, all contacts between the lowest two potential 
sheets occur as conic intersections, instead of tangencies, 
at points which are far from the potential minima. In par- 
ticular, both trigonal and pentagonal axes are locations 
of conical crossings for a in regions a and b respectively 
(i.e. when they do not correspond to minima). In partic- 
ular, for q on the % axis, the five electronic eigenvalues 
are given in Eq. (Q), where it can be readily verified that 
for a in region b the most negative one is indeed twofold 
degenerate. 

In region a, the six minima are all equidistant, defin- 
ing the simplest regular polytope in 5 dimensions (see 
Fig. [l]a). In this case, therefore, minimal closed paths 
join any of the 20 triplets of minima. It is straightfor- 
ward to verify that at the center of all such triplets there 
lies one of the trigonal axes, carrying a conical intersec- 
tion. If the degeneracy were restricted to the trigonal 
axes, however, the rich topology of the 5-D space would 
allow the triangular loop to squeeze continuously to a 
point avoiding the degenerate line: the associated Berry 
phase would then vanish. Instead, we checked that the 
two lowest sheets remain in contact through a bulky 3-D 
(1 radial + 2 tangential) region of distortions surrounding 
each trigonal axis. This guarantees the nontrivial topol- 
ogy of the loops, thus the possibility of nonzero Berry 
phase. Indeed a geometrical phase of it is associated to 
these triangular loops, as we computed explicitly by the 
discretized phase integral of Ref. JlSj] . Paths encircling 
two (or any even number) of such triangles (thus looping 
through 4, 6,... minima) have zero Berry phase, since 
the two phases cancel out. However, such paths, though 
energetically equivalent to the basic triangles (since they 
cross the same saddle points), are longer, therefore less 
relevant from a minimum-action point of view. We con- 
clude consequently that, for a in region a, the H ® h 
model must show the signature of a Berry-phase entan- 
glement. 

In region b, the minima are ten, each with 3 near- 
est neighbors and 6 second neighbors. The shortest 
closed paths through minima joins three points such as 
(l-»2->3->l) in Fig. |b. However, energetically, such 
loop is not the most convenient, since the segment join- 
ing two far neighbors (3— >1) must cross a barrier energet- 
ically 60% more expensive than that linking next neigh- 
bors (1— >2). Since the energy gaps between minima and 
saddle points grow as g 2 , eventually at strong coupling 
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only the "cheapest" paths affect the low-energy dynam- 
ics, and the relevant Berry phases should be calculated 
along such loops. Here, therefore, at large g, the low- 
energy paths are pentagons, such as (1^2^3^4^5^1) 
in Fig. [l]b. We computed the Berry phases for both kinds 
of paths, obtaining w and for the 3-points and 5-points 
loop respectively. This implies that the pentagonal loop 
encircles an even number (most likely 6) of degenerate 
regions (one of which around the pentagonal axis at the 
center of each 5-points loop), each carrying a phase fac- 
tor e l7r . We conclude that, in region b, although non- 
trivial Berry phases are present, they have no effect on 
the strong-coupling low-energy spectrum. Thus, in par- 
ticular, the GS symmetry should remain H in region a, 
while a nondegenerate A state must turn lower in region 
b at strong coupling. We stress that we have established 
the presence of nonzero Berry phases for all values of a, 
but also that, in region b, the effect of the geometrical 
phase is bypassed by energetically cheaper paths with 
null phase. 

This scenario is confirmed by numerical diagonaliza- 
tion (Lanczos method). ||J2(J I n Fig. |2| we plot the gap 
between the lowest H and A vibronic states, wherever 
Eh — Ea > 0, and where the GS is H. At weak cou- 
pling, as suggested by continuity, the GS is H. For g > 7 
and a in range b, the A state becomes the GS. We note 
however a little modulation in the boundaries of this re- 
gion, both g- and a-wisc. We observe, in particular, that 
the two special values a s and n — a s , far from mark- 
ing the closing of the H — A gap, show instead a rather 
sharp peak in the a direction. By drawing (in Fig. ||) the 
gap multiplied by g 2 , we evidence, along these ridges at 
a s and tt — a s , the g~ 2 large-g behavior of the H — A 
gap, characteristic of the motion in a flat trough of size 
~ g. Inside the region b, instead, the gap vanishes much 
more quickly, due to the tunnelling integral through the 
barriers between trigonal minima vanishing exponentially 
in g 2 . 

It is straightforward to extend the one-mode Hamil- 
4) to a more realistic case of many distortion 
2l| each characterized by its own frequency, cou- 



toman 
modes, 
pling and angle of mixing: 

H = } [-ffharm(^i) + H c _ v (gihuJi, ttj)] 



(8) 



We study in detail the two-modes case. Five free pa- 
rameters (u>i being taken as a global scale factor) ap- 
pear in the model. In order to carry out a significant 
study of the phase diagram, we limit ourselves to (i) two 
values only (1 and 5) of the ratio wa/^i, J22| and (ii) 
a-2 — n/2 = ai = a, assuming a principle of "maximum 
difference" between the modes. We take advantage of 
spectral invariance for individual sign change of each of 
the couplings gi — > — gi and for a — > —a, restricting 
to the < a < ir/2, g. L > sector. For convenience, 
we introduce polar variables g\ — <?cos7, <?2 = 3 sin 7 
(0 < 7 < 7r/2), and draw slices of the parameters space 
for fixed values of g, as a — 7 planes. 



The first interesting observation concerns the case of 
equal frequencies: even though Hamiltonian (^) is lin- 
ear in the coupling parameters the CG coefficients and 
the boson operators, the special case u>\ — L02 cannot be 
trivially reduced to a one-mode problem, by means of a 
suitable rotation mixing mode 1 and 2. This is a conse- 
quence of the linear independence of the coupling matri- 
ces V^ m '{a) for different values of a. 

We resort to exact diagonalization to treat the two- 
modes case. Due to the larger size of the matrices, we 
are limited to smaller couplings: we obtain a satisfacto- 
rily converged Eh — Ea gap up to g < 10 only. The 
calculations, for both uj^/wi = 1 and 5, show that for 
g < 7 the GS symmetry remains H for any a and 7 as 
in the one-mode case. Then, already at g — 8, an A 
(nondegenerate) GS makes its appearance in two local- 
ized regions of the a — 7 plane. Starting from g > 9, 
these separated regions assume essentially their asymp- 
totic strong-coupling shape (see Fig. ||). The first re- 
gion, located symmetrically across a = n/2, corresponds 
mainly to mode 1 with b-type (no-Berry) coupling: mode 
2 (Berry-phase entangled in this region) acts as a weak 
perturbation, incapable to change the GS symmetry for 
small enough 7. On the other side, the second region of 
A GS is located around a = 0: there, it is mode 2 who 
is responsible for the no-Berry phase coupling, mode 1 
acting as a weak perturbation, for 7 close enough to 7r/2. 
For u>2/u>i = 1 (not reported here), the two A-GS re- 
gions are, of course, equivalent. For = 5 (Fig. |^) 
instead, these two regions differ in size, in relation with 
the different relative energetics of mode 1 versus mode 2. 

In conclusion, we have illustrated the importance of 
the energetics of paths surrounding the points of degen- 
eracies of the two lowest BO potential sheets, for defining 
the effective role of the Berry phase. In all classical lin- 
ear JT models, the low-energy paths are affected by the 
geometrical phase in a way leading to a "boring" fixed 
ground-state symmetry. The H ® h model is special in 
being determined by an additional parameter, allowing to 
change the connectivity of the graph of low-energy paths 
through minima along with the regions of degeneracy of 
the two lowest sheets. Consequently, this new parame- 
ter leads continuously from a regular, Berry-phase entan- 
gled, region to a whole region where, although present, 
the Berry phase is totally ineffective in imposing its se- 
lection rules to the low-energy vibronic states, and to the 
GS in particular. 

Finally, for a system such as C^" our study implies 
that a detailed knowledge of not only the coupling pa- 
rameters gi, but also the characteristic angles on should 
be acquired for all modes in order to compute even such 
a basic property as the GS symmetry. 

We thank Arnout Ceulemans, Brian Judd, Fabrizia 
Negri, Erio Tosatti, and Lu Yu for useful discussions. 
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FIG. 1. The connectivity of the BO potential minima for 
regions a and b of angle a introduced in the text. All lines 
(solid and dashed) join nearest neighbor minima. 




FIG. 2. g 2 times the gap A = Eh — Ea (units of Tuo, 
logarithmic gray-scale) between the lowest non-G and non-A 
vibronic states as a function of g and a. In the positive re- 
gion, the GS is A, elsewhere it is H. This generates on the 
g — a plane a zero-temperature "phase diagram" . The basis 
is truncated to include up to 40 oscillator states. 




FIG. 3. The gap (units of %u)i) between the lowest non-G 
and non-^4 vibronic states as a function of a and 7 (defined 
in the text), for g — 10, UJ2/1-O1 = 5. In the positive region, 
the GS is A, elsewhere it is H. The basis includes up to 12 
oscillator states, enough to give a fairly converged value of 
A = E H -E A . 



4 



